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A Visualization of Null Geodesics for the Bonnor 
Massive Dipole 



Abstract 

In this work we simulate null geodesics for the Bonnor massive dipole metric 
by implementing a symbolic-numerical algorithm in Sage and Python. This 
program is also capable of visualizing in 3D, in principle, the geodesics for 
any given metric. Geodesics are launched from a common point, collectively 
forming a cone of light beams, simulating a solid-angle section of a point source 
in front of a massive object with a magnetic field. Parallel light beams also 
were considered, and their bending due to the curvature of the space-time was 
simulated. 


Resumen 

En este trabajo se simulan geodesicas nulas para la metrica de dipolo masivo 
de Bonnor, implementando un algoritmo simbolico-numerico en Sage y Python. 
Este programa es capaz, en principio, de visualizar en 3D las geodesicas de 
cualquier metrica. Estas geodesicas inician en un punto comiin, formando colec- 
tivamente un cono de rayos de luz, simulando una seccion de angulo solido de 
una fuente puntual frente a un objeto masivo con campo magnetico. Se consid- 
eran tambien rayos de luz paralelos, y se simulo su cambio de trayectoria debida 
a la curvatura del espacio-tiempo. 
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1 Introduction 

In General Relativity, Einstein field equations are a set of non-linear partial 
differential equations that describe the curvature of space-time due to the pres¬ 
ence of mass-energy-momentum. One of the exact solutions for the Einstein 
field equations in the vacuum was proposed in 1965 by W.B. Bonnor [T]; it de¬ 
scribes a massive spherical stationary object with a dipole magnetic field. This 
solution, expressed using the metric line element ds^ = g^ydx^dx^ Q, in which 
is the metric describing the manifold solution to these equations is: 

ds^ = + Zde^) “ ^ sin" dt^ (1) 

where 


P = — 2mr — b^cos^O (2a) 

Q = {r — m)^ — (6^ + im?)cos^0 (2b) 

Y = r‘^ — b^cos^O (2c) 

Z = — 2mr — b^ (2d) 


t^r^O^cj) are the standard spherical coordinates. The geometric system of units 
is used along this paper (c = G = 1). The parameters m and b are defined, 
according to Bonnor, as b being such that if r becomes large, the magnetic field 
is that of a magnetic dipole along 0 = 0 of moment mb7T~^ and the gravitational 
field becomes that of a mass 2m at the origin. [T] 

Trajectories on this metric must be found by using the geodesic equations, 
which can be thought as lines of minimum length (result from 5 f ds = 0): 

d?x'^ dx^ dx^ , , 

^Einstein summation convention applies. 


1 





subject to initial conditions that determine the kind of geodesic^ time-like 
{ds^ > 0, trajectories of particles with mass), light-like or null {ds^ = 0, tra¬ 
jectories of photons), or space-like {ds^ < 0, trajectories of hypothetical faster- 
than-light particles). A is an affine parameter, and are the components of 
the Christoffel symbol. 

In this paper null geodesics for the Bonnor metric are calculated in order 
to simulate the effect of the curvature of the space-time on a point source of 
electromagnetic radiation located at a certain distance from the origin, where 
the massive object is located. To achieve this purpose, a program in Python 
and Sage was developed to calculate analytically the Christoffel symbols from 
the metric and use them to build the geodesic equations. Once this is done, 
the equations are solved numerically using user-specified initial conditions, and 
these conditions are then varied to simulate a cone of light beams that represents 
a solid-angle section of the source. 

Other efforts have been made in this topic. For example. Miller and Frauen- 
diener (2010, [5]) developed Geodesic Viewer, which implements a selection of 
metrics which does not include the one studied here; and unlike this program, 
a local frame (tetrads) is required to specify the initial direction of the light 
beam. The justification for not using a local frame is that our scope is to study 
the collective behavior of geodesics, and high precision knowledge of the local 
directions is not required, as it would be for calculating gravitational leasing for 
such a point source. The Bonnor geodesics were also previously studied in 2D 
by Sanabria and Valenzuela [2], for charged test particles, as well as by Kovaf 
et al. [3], for neutral particles, from the point of view of orbit stability. 

Sage is an integrated environment consisting of several open-source packages 
(last count 89) with a common Python-based interface. It is capable of han¬ 
dling algebraic expressions (through Maxima and Sympy), as well as numerical 
calculations (through Scipy and Numpy), and many other tasks. 

Vpython (’’visual Python”) is a 3D modeling library for Python, that uses 
OpenGL for rendering graphics. Its main features are its simplicity, and real¬ 
time visualization capabilities. [6] 


2 Equations of motion for the Bonnor metric 


Expressions for the Bonnor geodesic equations were also calculated by means of 
the variational principle, to contrast them with the results generated with Sage. 
These expressions are: 
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^The sign of ds‘^ that determines kind of geodesic depends on the signature of the metric. 
In this paper, the signature of the metric is (1, —1, —1, —1) 
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Figure 1; Structure of the program 


Id fY^P^\ 
2 g3 j 


6»2 - — In 


dr 


y2 p2 

~W~ 


2y4 


3 r 


d 


d9 


e^^ln —r - 


Y2 


P d 
ZsiiPd~^ 


In 


'’’'+2ZW‘" 


sin^ 




p2 


( 5 ) 


where the dots mean derivative with respect to A, and e and I are integration 
constants. 


3 Our program 

3.1 Analytical parts of the program 

As sketched in the figure [TJ the program is formed by several modules. In 
principle, the program is capable of taking any metric written in any coordi¬ 
nate system and plotting geodesics with very little changes. The first module, 
equations . sage, computes the Christoffel symbols for the given metric analyt¬ 
ically using Sage, and stores the result as a Python script that can be later used 
to do calculations using only Python. This Python script, named by default 
geodesicEquations .py, contains the right-hand side of the Eq. [^defined as a 
function, as well as the metric. 

3.2 Numerical calculations 

The script geodesicEquations .py is used in the second module, geodesics .py 
to calculate and visualize the geodesics. The 4th-order Runge-Kutta method is 
used to solve the set of differential equations, and it is programmed in rk4.py. 
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The numeric calculations may also be handled using Numpy, a Python library 
with a variety of optimized functions written for this purpose; or a specific 
numerical method can also be integrated to the code this way. 

The initial conditions are given to this second module, and a function is called 
to vary them in order to form the light beam cone. Also, these initial conditions 
are varied in order to get a parallel set of beams. The given initial conditions are: 
the initial four-vector for the position, and the three spatial initial directions 
v [1] , v [2] , v [3] for the light beam; the other one is calculated by a function 
named calculate_vO() that uses the condition that ds‘^ = 0 for null geodesics. 
Since we are not using local frames, these spatial directions are really specified 
respect to a distant observer, assuming an asymptotically flat metric. Using a 
local frame and the condition of null geodesics, the spatial initial directions can 
be specified using only two angles, as done in [5]. 

Users can give input to the program by creating scripts that call the files 
equations. sage and geodesics .py. In this way, multiple metrics and initial 
conditions can be tested with no changes to the central code, and more func¬ 
tionality can be added in the future. 

The number of iterations and the step size can be specified, if not, default 
values are used. To avoid the effects of singularities on the behavior of geodesics, 
a stopping point based on distance to the origin is implemented and also can be 
specified in the main script. 

3.3 Visualization 

Vpython is used not only to plot the data, but also to manipulate the view 
in 3D in real time. Objects type curve () are used for the geodesics, and two 
small spheres mark the positions of both the source and the origin, where the 
massive object is located. A distinct arrow also points in the direction of -^z 
to indicate the direction of the dipole in the case of the Bonnor metric. While 
it is true that the code works for any coordinate system, a transformation to 
the Cartesian coordinate system is needed in order to represent positions in 
So far, transformations from spherical and cylindrical coordinates to Cartesian 
and vice-versa are already implemented. 

4 Results 

4.1 Test runs 

The program was tested first with Minkowski and Schwarzschild metrics. In 
the Fig. 1^ a light beam cone is drawn in a Minkowski space-time (in spherical 
coordinates), presented here to compare with the results obtained with other 
metrics. As it can be seen, each geodesic starts at the same point, marked by a 
small red sphere in the figure. As expected, geodesics are straight lines. Results 
for Schwarzchild metric are shown in Fig. |^. Fewer geodesics are displayed this 
time to better show how they are curved by the space-time. As expected, the 
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Figure 2: Test runs with point sources 



Figure 3: Test runs with parallel beams 
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Figure 4: Geodesics for the Bonnor metric with b = 0.001 
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light beams even cross each other when passing very close to the event horizon, 
and they are just slightly deviated if they pass farther away. 

Now, parallel beams are simulated for both metrics, and the results are 
shown in Fig. H In both cases, the direction of the beams is the negative x axis, 
and ah the geodesics remain in the x-y plane. The mass for the Schwarzschild 
metric is set to 0.5. 

The main scope for running these test metrics, where geodesics are known, 
was to detect possible errors and select more carefully the initial conditions for 
later runs. 

4.2 Point source with h <1 

Now, to the Bonnor metric. In Fig. nine geodesics have been simulated and 
viewed from different angles, to give a sense of spatial orientation. The source is 
locatec0at (10, 0, 0), and the direction of the central geodesic is radially inwards, 
varied then in —0.8 < v[2] < 0.8 and —0.3 < v[3] < 0.3, corresponding to the 
0 and (j) directions. 20 iterations with dlam = 0.1 (step size) have been carried 
out. The metric parameters are m = 1 and b = 0.001. 

The presence of the magnetic held makes the photons move in an helical 
motion around the ^ axis, that is, around the poles of the magnetic dipole, 
forming a jet-like structure. This behavior is consistent to the behavior of null 
geodesics for the Melvin metric, that is, locally, the magnetic held is dominant 
(see g]). 

4.3 The values of b and m 

Various values of b and m were tested. The behavior obtained in the preceding 
section corresponds to tested values 1 < 5 < 10“^^. With values of b in the 
interval 2 < 5 < 10 at least, the swirling seems to disappear. As predicted by 
Bonnor, with a value of m = 0 the space-time becomes hat and the result is 
similar to Fig. With a value of 6 = 0, also pointed out by Bonnor, the metric 
does not tend to the Schwarszchild one, but can be interpreted as a monopole of 
mass 2m together with higher mass multipoles dependent on m. Parallel beams 
of light on both metrics show the diherence (compare Fig. and Fig[^). 

4.4 Light beams with 6^5 

For the tested values of b = 2, 5,10 and m = 1 the results were similar. In 
Fig. §3, the same x-y view is shown, for comparison with Figs. and IDd. 
The magnetic held bends the light in such a way that a big opening is created 
around the object, and the light seems not to collide with it. Now, in Fig. 
variation in the z axis is tested to see the ehect on the distance for the magnetic 
dipole. As expected, the ehect decreases with distance, even in the z axis where 

^These position vectors are given in Cartesian coordinates; units are in the geometric 
system. 
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Figure 5: Geodesics for the Bonnor metric 



the dipole points to. It is important to notice that for Figs. and [^, the 
curves shown are all in the same plane. 

A point source located at (15, 0, 0) shows a different story (Fig. §i). Geodesics 
are no longer necessarily contained in a plane, and bend creating an ’’envelope” 
around the object that is less evident with the distance. 


5 Physical remarks 

Parallel beams of light simulate distant sources and therefore can be used for 
gravitational lensing, helping to diagnose which kind of object causes a partic¬ 
ular deformation on an image. Neutron stars can have very strong magnetic 
fields, up to the point of not being ignorable by General Relativity, but also 
have very rapid rotation rates, usually dominant. Light that behaves as shown 
in this paper, then, alerts of the presence of a massive object with a very strong 
magnetic field, but with negligible rotation. A point source, in the other hand, 
is also important to take into account, since several of these massive objects 
(black holes or neutron stars) also have companions, that is, nearby stars that 
act like point sources of radiation. 
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